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Abstract 

The  purpose  of  the  study  was  to  determine  the  relative  merits  of  the 
computer  codes  PF-CCMP  and  CAPS-2  for  calculating  radiation  fallout  pro¬ 
tection  factors  for  shelter  areas .  These  codes  were  produced  for  the 
Office  of  Civilian  Defense  (OCD)  by  the  Research  Triangle  Institute  (RTl) 
and  the  Architect-Engineering  firm  of  Praeger,  Karnaugh,  and  Waterbury, 
respectively. 

Protection  factors  for  various  detector  positions  within  five  building 
designs  of  varying  complexity  were  hand  calculated  using  the  Engineering 
Manual  method  as  outlined  in  the  most  recent  revision.  These  results  were 
then  compared  with  machine-calculated  protection  factors  using  the  PF-CCMP 
and  CAPS-2  programs.  As  a  result  of  the  comparisons,  some  program  errors 
were  found  in  PF-CCMP.  These  errors  were  corrected  by  the  authors  of  the 
code  and  the  cases  were  recalculated. 

The  PF-CCMP  code  was  found  to  calculate  protection  factors  in  the 
1-100  range  for  the  five  building  designs  to  within  +15$  of  the  hand-cal¬ 
culated  values .  The  protection  factors  calculated  with  the  CAPS-2  code 
tended  to  be  high  and  were  within  ~kk  to  +90$  of  the  hand-calculated  values . 

For  protection  factors  above  100,  the  PF-COMP  code  gave  results  that 
were  within  -4l  to  +36$  of  the  hand- calculated  values  and  tended  to  be 
conservative.  In  contrast,  CAPS-2  results  were  generally  not  conserva¬ 
tive,  with  the  percentage  error  ranging  from  -10  to  +58$.  Based  on  these 
calculations,  along  with  the  relative  ease  with  which  the  PF-COMP  code  may 
be  used  for  multiple-story  calculations,  the  PF-COMP  code  is  deemed  the 
better  code  and  the  use  of  CAPS-2  shouli  be  restricted  to  calculating 
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structures  with  low  protection  factors  (about  20  or  less)  and/or  of  only 

limited  complexity.  * 

I.  Introduction 

4 

The  Engineering  Manual  method  for  calculating  building  protection 
factors  involves  a  large  number  and  variety  of  design- type  equations  for 
which  much  of  the  shelter  configuration-dependent  data  are  available  in  the 
form  of  charts  and  graphs  of  limited  resolution.  The  synthesis  of  a  solu¬ 
tion  to  a  problem  of  average  complexity  is  a  tedius  and  oftentimes  confusing 
exercise,  and  a  single  calculation  for  a  given  problem  is  of  uncertain 
reliability  without  some  kind  of  corroboration  (which  is  usually  unavail¬ 
able). 

In  order  to  reduce  the  effort  required  to  analyze  structures  for 
shielding  effectiveness  against  fallout  radiation  and  the  possibility  of 

la 

error,  two  computer  codes  which  have  the  Engineering  Manual  method  as  a 

basis  were  produced  for  use  by  the  Office  of  Civil  Defense  (OCB).  There  j 

codes  are  FF-CCMP,^  produced  by  the  Research  Triangle  Institute  (RTI),  and 
p 

CAPS-2,  produced  by  the  Architect-Engineer  firm  of  Praeger,  Kavanaugh, 
and  Waterbury.  These  two  computer  programs  perform  the  numerical  equiva¬ 
lent  of  the  reading  of  charts  and  graphs  and  do  the  many  calculation 

sequences  which  would  be  involved  in  a  hand  calculation.  The  only  input  j 

required  is  the  configuration  data  for  the  shelter.  I 

Use  of  such  computer  programs  permit:  detailed  evaluations  of  many 
structures  for  a  small  cost  and  avoids  many  of  the  human  errors  that  would 

m 

be  prevalent  in  tedious  hand  calculation.  However,  the  usefulness  and 
reliability  of  the  computer-calculated  protection  factor  depends  largely 
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on  xxov  accurately  the  computer  code  represents  the  datQ,  equations,  and 
the  many  matters  of  judgement  which  comprise  the  Engineering  Manual  method. 
The  purpose  of  this  investigation  is  to  establish  how  well  the  computer- 
calculated  protection  factors  compare  with  hand- calculated  values  ‘using 
the  Engineering  Manual  method  for  a  vide  variety  of  shelter  configurations 
and,  on  the  basis  of  these  comparisons,  lo  achieve  an  evaluation  of  the 
PF-CCMP  and  CAPS -2  programs . 

The  procedure  followed  was  to  specify  a  set  ^f  hypothetical  building 
designs  which  are  representative  of  the  structures  that  would  be  encountered 
in  practice,  to  calculate  the  protection  factors  at  various  locations  with 
both  PF-CCMP  and  CAPS-2,  and  then  to  compare  the  results  with  the  hand- 

calculated  values  obtained  by  direct  application  of  the  latest  version  of 

3 

the  Engineering  Manual  method. 

In  the  following  sections  the  Engineering  Manual  method  find  the  two 
computer  codes  are  briefly  discussed.  The  results  of  the  calculations, 
both  by  hand  and  by  machine,  are  compared  in  graphical  and  tabular  form 
and  conclusions  are  drawn  from  these  comparisons.  In  the  appendix  a 
sample  calculation  for  one  detector  is  shown  in  order  to  demonstrate  the 
tediousness  of  the  calculation  and  how  the  Engineering  Manual  method  was 
applied. 


U.  The  Engineering  Manual  Method 

The  equations  and  basic  data  which  comprise  what  is  commonly  referred 
to  as  the  Engineering  Manual  method  were  developed  by  Eiserihauer  with  the 
help  of  L.  N.  FitaSimons  of  the  Office  of  Civil  Defense  from  the  basic 
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mot's,  by  Spencer. ^  The  details  of  the  method  are  given  in  a  number  of 

3-6 

publications . 

The  fundamental  approach  to  the  calculation  of  a  protection  factor 
(the  reciprocal  of  the  reduction  factor)  is  to  determine  the  reduction 
factor  associated  with  the  floors  above  and  below  the  detector  as  well  as 
the  detector  floor.  To  these  are  added  a  reduction  factor  due  to  the  roof 
and  mass  thickness  between  detector  and  roof,  known  as  the  "overhead  contri¬ 
bution"  to  the  reduction  factor.  The  floor  below  analysis  considers  the 
ground  direct  and  scattered  radiation.  The  detector  floor  includes  ground 
direct,  scattered,  and  skyshine,  while  the  floor  above  takes  into  account 
scattered  and  skyshine  contributions  . 

Equation  1  is  a  typical  equation  for  calculating  the  detector  floor 
reduction  factor  for  ground  direct,  scatter,  and  skyshine  contributions 
for  the  case  of  interior  partitions  and  no  mutual  shielding: 

°g  -  {C°d<“”Hd>  +  G.<V  '  P.  G.  -  S„,<Xe>]  Be<Xe>B>  +  Pa  Ga(“a> 

Be(0,B}  +  [Gs(%)  ♦  Gs(^)  -  Pa  Ga(a>a)3  SB(Xe)  E(e)  B,(Xe,H)} 

Bi(Xi)  ,  (1) 

where 

G^CtDijHj)  =  directional  response  for  ground  direct  contribution  to 
reduction  factor  based  on  solid  angle,  and  heigh., 

Hj,  above  contaminated  plane, 

Gs(®^)  =  directional  response  for  wall-scattered  ground  contribu¬ 
tion  to  reduction  factor  based  on  solid  angle  ti)4 , 

X 


G  (<u. )  =  directional  response  for  skyshine  ground  contribution  to 

reduction  factor  based  on  solid  angle  -j^, 

S^(Xg)  =  fraction  of  emergent  radiation  scattered  in  exterior  wall 

of  mass  thickness  X  , 

e* 

E(e)  =  shape  factor  for  wall-scattered  radiation  based  on  struc¬ 
ture  eccentricity  e, 

Be(Xg,H)  =  exterior  wall  barrier  reduction  factor, 

Pe  =  perimeter  ratio  of  apertures, 

I?i ( Xf )  =  interior  wall  barrier  reduction  factor  for  ground  contri¬ 
bution  based  on  interior  wall  mass  thickness  X^. 

The  overhead  contribution  considers  skyshine  through  the  roof,  as  well  as 
scattered  and  direct.  Methods  for  handling  apertures,  limited  fields  of 
contamination,  effects  of  interior  partitions,  and  the  use  of  fictitious 
buildings  to  determine  protection  factors  for  non-rectangular  buildings 
are  then  developed  to  complete  the  basic  framework  of  the  method. 


CAPS-2  and  PF-COMP  Computer  Codes 


The  computer  codes  CAPS-2  and  PF-CCMF  were  written  to  facilitate  the 
use  of  the  Engineering  Manual  method  by  eliminating  the  tedious  hand  calcu¬ 
lations  . 


CAPS-2 

The  CAPS-2  program  was  originally  developed  by  the  Architectural  and 
Engineering  firm  of  Praeger,  Kavanaugh,  and  Waterbury.^  Extensive  modifi¬ 
cations  were  made  by  Dirst  of  the  OCD  to  increase  the  flexibility  of  the 

o 

computations  and  the  input-output  routines. 


CAPS -2  is  written  in  FORTRAN  computer  language  and  is  based  on  the 
Engineering  Manual  method  with  one  important  exception  —  skyshine  contri¬ 
bution  for  the  detector  floor  and  the  floor  above  is  considered  to  be 
affected  by  mutual  shielding  buildings .  The  Engineering  Manual  method 
treats  this  contribution  as  the  same  with  or  without  mutual  shielding. 

The  code  will  handle  up  to  ten  detector  positions  per  floor  or  shelter 
area.  In  the  analysis  of  a  building,  the  program  logic  follows  that  of  a 
typical  engineering  manual  hand  calculation  and  includes  the  effects  of 
the  roof,  exterior  walls,  apertures,  interior  partitions,  floors,  mutual 
shielding,  height  above  contaminated  planes  and  building  geometry. 

Certain  restrictions  are  contained  in  the  program,  the  more  important 
ones  being:  (a)  window  sills  are  at  detector  level  (except  for  basement), 
which  is  3  ft  above  detector  floor,  only  two  contaminated  planes  may  be 
considered  for  each  side  of  a  structure,  and  (c)  the  exterior  wall  mass 
thickness  is  restricted  to  only  one  change  between  detector  floor  and  adja¬ 
cent  floors.  Preparation  of  the  input  data  is  relatively  simple  but  a  new 
set  must  be  prepared  for  each  detector  floor.  The  program  allows  for  a 
total  of  four  output  options  depending  on  the  desires  of  the  user.  The 
output  is  printed  in  a  manner  which  permits  an  analysis  of  those  areas  in 
which  the  effect  due  to  shielding  modifications  will  be  most  significant. 

PF-CCMP 

3.  0 

The  PF-CCMP  code  >  is  written  in  FORTRAN  computer  language  and  its 
use  is  restricted  to  the  large  computers  such  as  the  CDC  3^00.  The  pro¬ 
gram  is  also  based  on  the  Engineering  Manual  method  but  it  is  more  compre¬ 
hensive  than  CAPS- 2.  Roof  setbacks  are  included  and  allowances  are  made 
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I'or  three  contaminated  planes,  basement  areaways,  and  partial  basements. 
The  program  allows  for  different  interior  and  exterior  wall  mass  thick¬ 
nesses  for  each  floor.  In  addition,  sill  levels  are  not  restricted  to 
detector  level  (3  ft  fr an  floor),  but  only  one  change  is  permitted  above 
the  second  story.  Preparation  of  the  input  data  is  more  detailed  than 
CAPS-2,  but  it  need  be  done  only  once  for  a  building  since  PF-CCMP  will 
calculate  a  protection  factor  for  the  center  detector  position  of  each 
floor  in  addition  to  eight  other  locations  for  each  floor  (pre-set  by 
code).  The  output  is  essentially  the  same  as  CAPS-2. 


IV.  Building  Designs 

The  descriptions  of  the  five  hypothetical  buildings  (designated  as 
building no. i,building no. 2^tc.)  and  their  surroundings  are  presented  only 
in  the  detail  required  for  the  Engineering  Manual  hand  calculation.  Plane 
and  elevation  views,  dimensional  details,  construction  characteristics,  and 
detector  locations  are  given  for  each  building.  In  each  case  considered, 
the  detector  was  located  at  the  midpoint  and  3  ft  above  the  floor. 

Building  No.  1 

Building  No.  1  comprises  a  set  of  13  similar  buildings  (these  buildings 
are  designated  as  1A,  IB,  1C,  ...,  1M)  for  which  only  minor  design  differ¬ 
ences  exist.  The  data  generated  for  this  case  allowed  more  detailed  compar¬ 
isons  to  be  effected  wherein  the  influence  of  one  portion  of  the  overall 
calculation  could  be  determined.  The  plan  and  elevation  views  are  shown 
in  Fig.  1,  and  the  design  specifications  are  given  in  Table  1.  A  protec¬ 
tion  factor  was  computed  for  the  center  position  of  the  detector  in  the 
basement  and  in  the  second  floor. 
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Fig.  1.  Building  No.  1  Plan  and  Elevation  Views. 
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Building  Mo.  2 

The  plan  and  elevation  views  of  building  no.  2  are  shown  in  Fig.  2. 
The  shading  in  the  plan  view  indicates  differed  roof  elevations  and  the 
numbers  in  the  center  of  a  roof  designate  the  total  height  of  that  sec¬ 
tion  above  ground  level. +  Design  specifications  are  presented  in  Table  2. 
A  protection  factor  was  computed  for  the  first  and  fifth  floors  with  the 
detectors  centrally  located. 

Building  No.  3 

Building  no.  3  is  typical  of  many  commercial  buildings  such  as  small 
factories  or  hospitals .  The  plan,  east-wall  elevation,  and  south-wall 
elevation  views  are  shown  in  Fig.  3*  Partition  layouts  for  all  four 
floors  are  shown  in  Fig.  4.  Design  specifications  are  presented  in 
Table  3*  Protection  factors  were  computed  for  the  first  and  fourth  floors 
with  the  detectors  centrally  located. 

Building  No.  4 

Building  no.  4,  together  with  its  surroundings,  is  typical  for  a 
large  apartment  or  office  building  in  a  metropolitan  area.  The  plan  and 
elevation  views  are  shown  in  Fig.  5,  and  the  partition  layout  for  all 
floors  is  shown  in  Fig.  6.  The  design  specifications  are  given  in 
Table  4.  Protection  factors  are  computed  for  the  first  and  third  floors 
only. 

Building  Ho.  5 

Building  no.  5  could  be  a  part  of  a  hospital  or  school  complex  in  a 
suburban  area.  The  structure  is  not  synsnetric  as  was  the  case  for 

+This  comment  applies  to  all  plan  views . 


Dimensional  Data: 


Total  height,  60  ft 
Basement  height,  12  ft 
First  floor  height,  10  ft 
Upper  floor  height,  10  ft 

Aperture  heights,  k  ft  (bottom  to  top  of  aperture) 
Sill  height,  3  ft 

Construction  Specifications: 

Basement  wall  mass  thickness,  50  pst 
First  floor  wall  mass  thickness,  50  psf 
Upper  floor  wall  mass  thickness,  50  psf 
Basement  floor  mass  thickness,  35  psf 
First  floor  mass  thickness,  4o  psf 
Upper  floor  mass  thickness,  35  psf 
Roof  mass  thickness,  45  psf 

No  interior  partitions 

Apertures : 

Basement,  Of> 

First  floor,  3Cff> 

Upper  floor,  30% 
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Fig*  3*  Building  No.  3  Plan  and  Elevation  Views 


<*♦  AND  2nd  FLOORS 


4.  Partition  Layouts  c~  Floors  1  to  4  in  Building  No.  3 
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Table  3.  Specifications  for  Building  No.  3 

Dimensional  Data: 

Total  height.  60  ft 
First  floor  height,  10  ft 
Upper  floor  height,  10  ft 

Aperture  height,  5  ft  (bottom  to  top  of  aperture) 
Sill  height,  3  ft 

Construction  Specifications: 

First  floor  wall  mass  thickness,  60  psf 
Upper  floor  wall  mass  thickness,  60  psf 
First  floor  mass  thickness,  30  psf 
Upper  floor  mass  thickness,  30  psf 
Main  roof  mass  thickness,  40  psf 
Setback  roof  mass  thickness,  40  psf 
Interior  partition  mass  thickness,  25  psf 

No  interior  partitions  on  4tn  and  5th  floors 

Apertures: 

First  floor  and  second  floor: 

North  and  south  walls,  35$ 

East  and  west  walls,  25$ 

Third  and  fourth  floors: 

North  and  south  walls,  35$ 

East  and  west  walls,  42$ 

Fifth  and  sixth  floors: 

North  and  south  walls,  35$ 

East  and  west  walls,  35$ 
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Fig.  5.  Building  No.  4  Plan  and  Elevation  Views. 


Dimensional  Data 


Total  height  of  building,  484  ft 
Namber  of  stories,  40  with  basement 
Basement  height,  l4  ft 
First  floor  height,  l6  ft 
Upper  floor  height,  12  ft 

Aperture  height,  1st  floor,  12  ft  (bottom  to  top 
of  aperture) 

Upper  floor  aperture  height,  8  ft 
Sill  height,  all  floors,  3  ft 

Construction  Specifications: 

Basement  wall  mass  thickness,  75  psf 
First  floor  wall  mass  thickness,  75  psf 
Upper  floor  wall  mass  thickness,  75  psf 
All  interior  partition  mass  thickness,  25  psf 
Basement  floor  mass  thickness,  80  psf 
First  floor  mass  thickness,  80  psf 
Upper  floor  mass  thickness,  35  psf 
Roof  mass  thickness,  70  psf 

Apertures: 

Basement,  0$ 

First  floor,  60 $ 

Upper  floors,  25$ 
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buildings  no.  1  through  no.  4,  and  its  surroundings  include  a  small  lake 
and  two  adjacent  buildings  of  comparable  size.  The  plan  view  of  the 
building  and  its  surroundings  is  shown  in  Fig.  7.  The  north-wall,  east- 
wall,  south- wall,  and  west-wall  elevations  are  shown  in  Fig.  8,  along  with 
the  details  of  an  areaway  adjacent  to  building  no.  5.  Partition  layouts 
for  all  five  floors  are  shown  in  Fig.  9 •  The  plan  view  of  the  basement 
is  shown  in  Fig.  10.  Design  specifications  are  given  in  Table  5.  Protec¬ 
tion  factors  were  calculated  for  the  partial  basement,  first,  and  fifth 
floors . 


Calculations!  Procedures  and  Comparison  of  Results 


Protection  factors  for  the  five  hypothetical  buildings  described  in 
the  previous  section  were  calculated  by  hand  using  equations  and  data 
found  in  TR-20.  These  calculations  were  performed  with  great  care  and 
the  results  should  represent  the  application  of  the  Engineering  Manual 
method  to  an  accuracy  within  the  limitations  imposed  by  reading  the  charts 
and  the  use  of  good  judgement  in  the  areas  of  uncertainty  in  application 
of  the  method.  The  aim  was  tu  compare  computer- calculated  protection  fac¬ 
tors  with  these  "standard"  values.  In  the  interpretation  of  the  results, 
the  basic  assumption  was  that  any  disparities  were  due  to  an  imperfect 
representation  of  the  Engineering  Manual  method  by  the  computer  codes. 


As  an  example  of  the  tedious  calculations  required  in  applying  the 
Engineering  Manual  method,  the  very  extensive  and  complex  details  of  the 
hand  calculation  of  a  protection  factor  for  one  detector  position  in  a 
multiple-story  building  (building  no.  5)  with  complicated  geometry  and 
mutual  shielding  are  given  in  an  appendix  to  this  report. 
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Fig.  8.  Building  No.  5  Elevation  Views  and  Areaway  Detail 


Dimensional  Data 


Basement  height,  12  ft 

First  floor  height,  10  ft 

Upper  floor  height,  10  ft 

Sill  heights  in  basement  areaway,  6  ft 

Aperture  height  in  basement  areaway,  4  ft  (bottom  to  top  of  aperture ) 
Upper  floor  aperture  height,  5  ft 
Upper  floor  sill  height,  3  ft 

Note:  On  the  south  wall,  the  first  contaminated  plane  is  4  ft 

lower  than  the  other  three  walls  over  the  entire  azimuthal 
sector  seen  by  the  south  wall. 

Construction  Specifications: 

Basement  wall  mass  thickness,  60  psf 

First  floor  wall  mass  thickness: 

North  wall,  50  psf 
East  wall,  50  psf 
South  wall,  50  psf 
West  wall,  60  psf 

Upper  floor  wall  mass  thickness  (2nd  and  3rd  floors): 

North  wall,  50  psf 
East  wall,  50  psf 
South  wall,  50  psf 
West  wall,  60  psf 

Upper  floor  wall  mass  thickness  (4th  and  5th  floors): 

North  wall,  45  psf 
East  wall,  45  psf 
South  wail,  45  psf 
West  wall,  45  psf 

All  interior  partition  mass  thickness,  25  psf 

Basement  floor  mass  thickness,  50  psf 

First  floor  mass  thickness,  4o  psf 

Upper  floor  mass  thickness,  35  psf 

Upper  floor  mass  thickness  (if  change),  30  psf 

Story  of  change,  4th 

Setback  roof  mass  thickness,  50  psf 

Main  roof  mass  thickness,  50  psf 


. in,, . j  o  it  m&sxmm .  — i - 
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Table  5-  (cont.) 


Apertures : 

Percentage  apertures  in  basement  areaway,  2C% 

Total  aperture  width  of  30  ft  in  wall  of  length,  50  ft 

Apertures  on  first  floor; 

North  wall,  25/6 
East  wall,  35/6 
South  wall,  3056 
West  wall,  4oj6 

Apertures  on  upper  floors  (no  change); 

North  wall,  25/6 
East  wall,  35% 

South  wall,  30$ 

West  wall,  hCFf, 

Apertures  on  upper  floors  (if  change): 

North  wall,  2556 

East  wall,  20%  (story  of  change,  kth) 

South  wall,  1556 
West  wall,  2056 


The  PF-COMP  and  CAPS -2  computer  codes  were  used  to  calculate  protec¬ 
tion  factors  for  the  same  five  buildings  (as  described  in  the  previous 
section)  using  exactly  the  same  input  data  whenever  possible.  The  CAPS-2 
computer  calculations  were  run  on  the  IS  I  l6o4  at  the  Oak  Ridge  National 
Laboratory  and  the  PF-CCMP  calculations  were  run  at  the  National  Civil 
Defense  Computer  Facility,  Washington,  D.  C.  The  operation  of  both  codes 
seemed  routine;  however,  a  detailed  comparison  of  the  PF-CCMP  results  with 
the  standard  (hand  calculated)  protection  factors  suggested  that  the 
PF-CCMP  code  was  in  error.  The  authors  of  the  code  were  informed  of  the 
discrepancies  and  they  were  able  to  locate  and  correct  the  errors,  which 

Q 

were  in  floor-above  and  floor-below  contributions .  The  complete  set  >1 


problems  was  recalculated  with  the  corrected  code  and  better  agreement 
was  obtained. 


Comparisons  of  the  calculated  protection  factors  are  presented  in 
graphical  form  in  Figs.  11  and  12.  Fig.  11  is  a  plot  of  the  protection 
factors  calculated  with  the  PF-CCMP  code  versus  the  corresponding  hand- 
calculated  values .  A  computer-calculated  protection  factor  would  exactly 
equal  the  hand-calculated  value  if  the  plotted  point  would  lie  on  the 
line  drawn  through  the  origin  at  ^5°.  A  point  located  below  the  line 
indicates  an  underestimate  of  the  protection  factor  by  the  computer  code 
with  respect  to  the  hand-calculated  /ilue;  a  point  tbove  the  line  indi¬ 
cates  an  overestimate.  Note  that  the  PF-COMP  values  generally  are  in  good 
agreement  with  the  hand  calculations  for  buildings  having  protection  fac¬ 
tors  up  to  about  100  and  generally  lie  well  below  the  ^5°  line  for  larger 
protection  factors . 

Fig.  12  is  a  plot  of  the  protection  factors  calculated  with  CAPS-2 
versus  hand- calculated  values.  The  CAPS-2  data  points  in  general  lie 
above  the  45°  line,  indicating  that  CAPS-2  generally  overestimates  the 
protection  factor  in  comparison  with  the  hand  calculations .  Note  that  the 
CAPS-2  data  points  form  a  rather  scattered  pattern  as  Caspar ed  with  the 
PF-CCMP  calculations  shown  in  Fig.  11. 

The  protection  factors  for  all  detector  locations  are  also  presented 
in  Table  6.  No  attempt  was  made  to  estimate  error  or  confidence  limits 
for  the  individual  protection  factors  because  of  the  nature  of  the  calcu¬ 
lations;  the  percentage  deviation  shown  is  relative  to  the  hand  calcula- 
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Fig.  11.  Comparison  of  Protection  Factors  Calculated  with  the 
■OF-CCMP  Code  and  with  the  Engineering  Manual  Method. 
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Fig.  12.  Comparison  of  Protection  Factors  Calculated  with  the 
CAPS -2  Code  and  with  the  Engineering  Manual  Method. 
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Calculated  Protection  Factors  Deviation  from  Haiti  Calc. 


Figs.  13  and  l4  are  riots  of  the  ''floor-above"  and  "floor-below" 
contributions  as  calcu laced  with  the  original  and  corrected  versions  of 
PF-C(3»!P,  respectively,  versus  the  corresponding  values  calculated  by  hand. 
These  plots  reveal  the  nature  of  the  errors  contained  ir.  the  original 
PF-CCMP  code;  the  original  version  of  PF-COMP  seriously  overestimated  the 
reduction  factors  in  most  cases  for  the  "floor-above"  and  "floor-below" 
contributions .  The  excellent  correlation  of  the  modified  PF-COMP  results 
with  the  hand-calculated  ones  generally  confirms  that  the  errors  in  the 
original  PF-COMP  code  have  been  resolved. 

VI.  Conclusions  and  Recommendations 
An  inspection  of  Fig.  11  and  Table  6  shows  that  the  protection  fac- 

Q 

tors  <  100  calculated  with  the  PF-CCMP  code0  are  in  good  agreement  (-12  to 

+17$)  with  the  hand-calculated  values  and  generally  tend  to  be  a  little 

conservative,  a  desirable  characteristic  for  a  calculation  of  this  kind. 

2 

In  contrast,  the  protection  factors  calculated  with  CAPS-2  code  tend  to 
be  high  and  do  not  agree  as  well  with  the  hand-calculated  values  (-44  to 
+90$).  The  deviations  from  the  hand-calculated  values  became  larger 
(-41  to  +36$  for  PF-CCMP  and  -10  to  +58$  for  CAPS-2)  for  protection  fac¬ 
tors  >  ,100.  For  this  range  only  PF-COMP  should  be  used  since  it  tends 
to  gi\e  conservative  values  in  contrast  to  CAPS-2. 

As  a  result  of  this  study,  along  with  the  relative  ease  wi„b  which 
PF-CCMP  may  be  used  for  multiple-story  calculations,  the  PF-COMP  code  is 
deemed  to  be  the  better  code  and  the  use  of  CAPS-2  should  be  restricted 
to  calculating  structures  with  low  protection  factors  (about  20  or  less), 
and/or  of  only  limited  complexity. 
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Fig.  13*  Comparison  of  the  "Floor -Be low"  Contributions  Calculated 
with  FF-CGMP  Codes  (Original  and  Revised)  and  with  the  Engineering  Manual 
Method. 
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Fig.  1^.  Comparison  of  the  "Floor- Above"  Contributions  Calculated 
with  PF-CCMP  Codes  (Original  and  Revised)  and  with  the  Engineering  Manual 
Method . 
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APPENDIX 


Sample  Calculation  for 
Building  No.  $  First  Floor 
Detector  -  Engineering 


Manual  Method 


Table  Al.  Equations  for  Calculating  Various  Contributions  to 
Reduction  Factor  by  Engineering  Manual  Method+ 
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floor;  no  mutual  shielding. 


Table  Al.  (cent. 


Since  many  of  the  equations  vised  in  the  engineering  method  are  used  repeatedly  in  the  solution  to 
this  problem,  this  table  was  prepared  for  ready  reference.  When  these  equations  are  used  in  the 
solutions,  they  are  referred  to  by  the  number  in  the  table. 


Ground  Contribution  Through  North  Wall  (see  Fig.  Al)+ 


Parameters  and  functions: 
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1.  Contribution  through  shielded  sector,  A  =  0.147  (North  Wall) 

Z1 

(a)  Detector  floor  (Eqs.  8  +  10  +  11  +  12,  Table  Al). 


+The  charts  referred  to  throughout  this  calculation  are  the  standard 
charts  used  in  the  Engineering  Manual  Method. 
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Cg  =  |(0.30-0. 175  )(0. 1*15  )(o.304)  +  [0.076+0.167-0.506(0.122)] 

x  (o.58)(l.i*)(o.l5)  +  [0.0U5-0. 506(0. 033)](0. U5X0.30I*) 

+  0. 506(0. 033)(l.0)j  (0.5i*)(0.ll*5)  =  0.00U3 
(b)  Floor  above  (Fqs .  14+15+  1 6,  Table  Al) 

C  =  {(0.074-0. 0l*5)(0.!*15)(0.30!*)(0. 75)  +  (0. 074-0. 0l*5)(0.25)(l.0) 

+  (0.31-0. l67)(0.56)(l.4)(0.15)(0.75)}  (0.54)(0.125)(0.l475) 

=  0.00023 

2.  Contribution  through  unshielded  sector,  A  -  A  =  0.0573 

z2  zi 

(a)  Detector  floor  (Eqs.  7  +  9  +  11  +  12,  Table  Al) 

Cg  =  -[[0.3+0.01*5-0.506(0.033)]  (0.1*15)(0.30l*)  +  0. 506(0. 033)(l.O) 

+  [0.076+0.167-0.506(0.122)](0.58)(1.4)(0.30!0}  (o.54)(o.0573) 
=  0.00313 

(b)  Flocr  above  (Eqs.  13  +  15  +  16,  Table  Al) 

Cg  =  -[(0.07  4-0. 0l*5)(0.4l5  )(0.304)  (0.75)  +  ( 0.C74-0.045)  (1.0)  (0.25) 
+  (0. 31-0. 167) (0. 585 )(l.4)(o.304)( 0.75)}  (0.54)(0. 125) (0.0588) 
=  0.00014 

Total  contribution  through  north  wall: 


0.00463  +  0.00023  +  0.00318  +  0.00014  =  0.00818. 


1*3 


Ground  Contribution  Through  East  Wall  (see  Fig.  A2) 
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1.  Contribution  through  shielded  sector,  A  =  0.135 

Z1 

(a)  Detector  floor  (Eqs .  8  +  10  +  11  +  12,  Table  Al) 

Cg  =  {(0.3-0. n)(o. 1*15X0.30**)  +  [0.076+0.167-0.7(0.122)]  (0.58) 

x  (1.10(0.199)  +  [0.01*5-0.7(0.033)]  (0A15X0.30I*)  +  0.7(0.033) 
x  (1.0)}  (0.5l*)(0.135)  =  0.00328 
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A2.  Geometry  for  Calculating  Ground  Contribution  Through  Sast 


(b)  Floor  above  (Eqs .  l4  +  15  +  1 6,  Table  Al) 

Cg  =  {(0. 074-0. 045)(0.4l5)(0.304)(0. 65)  +  (0*074-0. o45)(l.o)(0. 35) 
+  (0.31-0. 167 )(0. 585 )(l.4)(0.199)(o. 65 )j-  (0.125) (0.5^) (0.135) 
=  0.00025 

2.  Contribution  through  unshielded  sector,  A  -  A  =  0,l6 

Z2  Z1 

(a)  Detector  floor  (Eqs.  7+9+11+12,  Table  Al) 

cg  =  {[0.3+0.0^5-0.7(0.033)]  (0.4l5)(0.?04)  +  0.7(0.033)(1.0) 

+  [0.076+0.167-0.7(0.122)]  (0.58)(1. 4)(0.30M}  (0.5U)(o.16) 

=  O.OO887 

(b)  Floor  above  (Eqs.  13  +  15  +  l6,  Table  Al) 

Cg  =  {(o. 074-0. G+5)(o.Ui5)(o.30U)(o. 65)  +  (0. 074-0. 045)(1.0)(0. 35) 
+  <c  1-0. 167)(0.58)(1.4)(0. 3^0(0.65))  (0.125) (0.54) (0.16) 

=  O.OJ038 

Total  contribution  through  east  wall: 

0.005-48  +  0.00025  +  0.00817  =  0.01498 
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Grcur.d  Contribution  Through  South  Wall  (see  Fig.  A3) 
Parameters  and  functions: 
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0.14 

0.35 

0.35  j 

% 

75 
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10 

0.75 

0.20 

0.80  0.43 

0.235 

1 

f 

t 

a«s 

10 

6o 
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0.1667 

O.366 

oos  -  O.26/2  = 

0.13 

1 

1 

1 

^s 

10 

6o 

7 

0.1667 

0.233 

uig  =  0.40/2  = 

0.20 

i 
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a.  Chart  5* 

b .  Chart  6 . 


Be(60',7')  =  0.20 

(Chart  2) 

Be(50,7’)  =  0.26 

(Chart  2) 

Be(0.7’)  =  0.88 

(Chart  2) 

Be(60,ll')  =  0.18 

(Chart  2) 

Be(0,ll’)  =  0.8 

(Chart  2) 

s  (50)  «  0.58 

U) 

(Chart  7) 

3  (60)  =  0.63 

U) 

(Chart  7) 

B^(35)  =  0.125 

(Chart  l) 

E(.75)  =  1.4 

A  =  0.30] 
y  1 1st  Floor 

pa  =  0.6  J 

(Chart  8) 

i 
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A  =  0.20,  basement 

P 

B  (0.13,50)  = 

(JJS 

0.005 

(Chart  9) 

B'  (0.20,50)  = 

OB 

0.0128 

(Chart  9) 

Bi(25)  =  0.546 

(Chart  l) 

Bf(40)  =0.11 

(Chart  1) 

1.  Contribution  through  A  =  0.204 

Z1 

(a)  Detector  floor  (Eqs .  7  +  11  +  9  +  12,  Table  Al) 

=  {[0.2+0. 045-0.6(0. 033)]  (0.4i5)(o.26)  +  0.6(0. 033)(0.88) 

+  [0.C76+0.l67-0.6(0.122)]  (o.58)(l.4)(o.26)j  (0.546)(0.204) 
=  0.00854 

(b)  Floor  above  (Eqs.  13  +15  +  l6,  Table  Al) 

Cg  =  {(0. 074-0. 045)(0.4l5)(0.26)(0.7)  +  (0.074-C.045)(0.88)(0.3) 
+  (0.31-0.l67)(0.585)(l.4)(0.26)(0.7)j-  (0.546)(0.126)(o.2o4) 
=  0.00043 

2.  Contribution  throyerh  A  =  0.140 

“  Z2 

(a)  Areaway  sources 
(l)  Direct 

Cg  =  CGd^11’)  "  GdW,’11’)^1  '  Su>(Xe):i  Be(6°^11,) 
x  Bf(Xf)A  =  (0.47-0. 35)(0.37)(0.l8)(0.1l)(0.l40) 


0.00012 


. . . . . .  wwKMf  oy.jw’WMaagcAirtiiwea^fsgigHiit^  jyK^uHEWW«BSBWM<WUg!WBg 
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(2)  Scatter 


Cg  =  [GsK}  -  °a^3  W  E(e>  B..<VXe>  Bf<Xf>  % 

«  [0.245-0. 076](0.63)(1.U)(0.005)(0.11)(0.140) 

=  0.000012 

(b)  For  infinite  field  sources 
(l)  Direct 

Cg  -  {[Gd(u^,7')  -  Gd(ut,7’)][l  -  SjXe)]  Be(60,7’)  Bf(Xf) 
x  (1  -  Ap)  +  [G4(V7')  -  GdUL/7’)]  Be(0,7')  Ap} 


x  Bf(Xf )  kZr 


C6  ’  {(G« 43-0. 20) (0.37 )( 0.2) (0.7)  +  (0.43-0.20)(0.88)(0.3)} 


x  (0.126)(0.l4o)=  0.00.112 


(2)  Scatter 


ce  -  tGsk>  -  G>L)J  Vx«>  E<e>  VXP  \ 

=  (0. 215-0. 076)(0.63)(1.4)(0.0128K0.11)(0.140) 

=  0.000024 

Total  contribution  through  south  wall: 

0.00854  +  0.00043  +  0.00012  +  0.000012  +  0.00112  +  0.000024  »  0.01024 
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Ground  Contribution  Through  West  Wall  (see  Fig-  A4) 
Parameters  and  functions: 
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a. 

Chart  6. 

b. 

Chart 

5- 

Be(60,3’) 

-  0.238 

(Chart  2) 

Be(0,3*) 

=  1.0 

(Chart  2) 

s  (60)  ■ 

«u 

0.63 

(Chart  7) 

B'(35)  - 

0.125 

(Chart  l) 

B«* 

(( 0.46, 60)  ■ 

0.109 

(Chart  9) 

*(0.75)  - 

1  1.4 

(Chart  8) 

A  *  0.40 
P 

P  *  0.8 
a 

Bi(25)  =  0.54  (Chart  l) 

A_  =  0.295 
"1 

Detector  floor  (Eqs  .8  +  10  +  11  +  12,  Teble  Al) 


1. 


C  =  {(0.3-0. 195)(0.37)(0. 238)  +  [0.076+0.167-0.8(0.122)]  (0.63) 

X  (1. *0(0.109)  +  [0.0^5-0.8(0.033)]  (c.37)(0.238) 

+  0.8(0. G3^)(l.0)|  (o.5J0(o.295)  =  O.OO816 

2.  Floor  above  (Eqs .  1*4-  +  15  +  16,  Table  Al) 

Cg  =  |(0.07^-0.045)(0.37)(0.238)(0.6)  +  (0. 07^-0. 0^5)(i.0)(0J+) 

+  (0.31-0.167)(0.63)(1.4)(0.109)(0.6)|  (o.5*0(o.295)(o.l25) 

=  0 . 000*4-3 

Total  contribution  through  west  wall: 

0.00816  +  0.00043  =  0.00859 

Overhead  Contribution  (see  Fig.  A5a) 

The  overhead  contribution  consists  of  that  from  the  30-ft  high  set¬ 
back  and  the  50-ft  h.'gh  main  roof.  These  contributions  may  be  found  ho¬ 
using  the  fictitious  building  concept. 

The  setback  contribution,  which  is  the  contribution  from  the  shaded 
area  CQ^  (see  Figs.  A5b  and  c)  in  the  actual  building,  is  obtained  by 
subtracting  out  the  contribution  of  the  unshaded  area  from  that  obtained 
for  the  entire  roof  of  the  fictitious  building. 

X  =  120  psf 
o 

Bi(Xi)  =  Bj_(25)  =  0.^3  (Chart  l) 

«0  *  “>  3) 

«;=■»(  I5  ,  )  =  0.19  (Chart  3) 


ORNL-DWG  68-7133 
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Fig.  A5.  Geometry  for  Calculating  Overhead  Contribution  (Setback) . 
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Co(u,o,120)  =  0.0106  (Chart  4) 

C' (^,120)  =  0.0056  (Chart  4) 

Ccl  =  (0.6)  [CoUo,120)  .  Co(«b^120)]  BJ(X’) 

-  (0.6)  f 0.0106 -0.0058]  (0.43)  =.-  0.00124 

The  contribution  from  the  main  roof  is  more  difficult  to  calculate. 
Referring  to  Fig.  A5a,  partition  A  lies  over  the  detector  and  extends 
from  the  fourth  floor  to  the  roof;  hence,  it  may  be  neglected.  Partition 
B  extends  from  the  first  floor  to  the  fourth  floor  and  has  e  little  more 
effect,  but  it  still  is  too  small  to  warrant  the  added  complexity  of 
including  it.  The  only  partition  which  affects  the  main  roof  contribu¬ 
tion  significantly  is  partition  C,  which  shields  the  contribution  from 
the  shaded  area  Cq2  (see  Figs.  A6a  and  b).  This  contribution  is  obtained 
by  subtracting  the  unshaded  area  from  the  total  for  the  fictitious 
t”ilding  and  proportioning  on  an  area  basis; 

XQ  ■  180  psf 

B^(25)  =  0.43  (Chart  l) 

%  -  ■  <  $5  -  ^  >  -  °-30  !chart 

«£  -  «  (  .  ^22  )  =  0.19  (Chart  3) 

C  (0. 30,180)  =  0.0021  (Chart  4) 

C^(0. 19,180)  =  0.0016  (Chart  4) 

co2  ■  (C. 0021-0. 00l6)(0. 43)  =  0.00014 

The  remaining  portion  from  the  main  roof  (see  Figs.  A6c  and  d)  is 
s  traightf orward; 


3 


5  6 

=  180 

o 

B,*(XJ)  =  B'(25)  *  0.43  (Chart  l) 

A  A  1 

“o  *  w  (  7§  »  "7^  )  -  °-19  (Chart  3) 

CQ(0. 19,180)  =  0-0016  (Chart  4) 

co3  -  Co(0.19,180)  =  (0- 0016 )( 0.667)  =  0.00106 

The  total  overhead  contribution  to  the  reduction  factor  from  both  the 
main  roof  and  the  setback  is: 

C  =C.  +  C,+C,s  0.00124  +  0.00014  +  0.00106  =  0.00244 
O  ol  02  03 

Total  Protection  Factor  for  the  First  Floor 

R.F.  «=  North  Wall  Contribution  +  East  Wall  Contribution  +  South  Wail 
Contribution  +  West  Wall  Contribution  +  Overhead  Contribution 
=  0.00818  +  0.01498  +  0.01024  +  O.OO859  +  0-00244  =  0.0444 
or 

,  P.P.  «  23 
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The  relative  merits  of  two  coape ter  codes,  FF-CCKP  and  CAPS-2,  for  calculating  radiatioi 
fallout  protection  factors  for  shelter  areas  vere  investigated  by  comparing  the  code 
results  vith  those  froa  hand  calculations  based  on  the  Engineering  Manual  method.  Five 
building  types  vere  considered.  For  protection  factors  in  the  range  of  1  to  100,  the 
PF-CCKP  code  was  found  to  yield  values  that  vere  vithin'+15J&  of  the  hand-calculated 
values ,  vhile  the  CAPS-2  code  gave  results  that  v<  -e  vithin  -4b  to  +90J6  o'"  the  hand- 
calculated  values.  For  protection  fjctors  greater  than  100,  the  PF-CCMP  code  results 
vere  vithin  -4l  to  +3&fc  and  the  CAPS-2  code  results  vere  vithin  -10  to  +5 &f>  of  the  hand- 
calculated  values.  Bass.,  on  tfc^se  calculations,  along  vith  the  relative  ease  vith  vhicl 
tl.a  PF-CCMP  code  may  be  used  for  multiple-story  buildings,  the  PF-CCMP  code  is  deemed 
the  better  code  and  it  is  recommended  that  the  CAPS-2  code  be  restricted  to  structures 
vith  lov  protection  factors  (about  20  cr  less )  and/or  limited  complexity. 
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